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Abstract: Motivated by the findings of the OPERA experiment, we discuss the hypoth- 
esis that neutrino propagation does not obey Einstein special relativity. Under a minimal 

set of modifications of the standard model Lagrangian, wc consider the implications of non 
standard neutrino propagation on the description of neutrino interactions and, specifically, 
on the pion decay processes. We show that all the different dispersion relations which 
have been proposed so far to explain OPERA results, imply huge departures from the 
standard expectations. The decay channel 7r+ — )• e'^u^ becomes significantly larger than in 
the standard scenario, and may even dominate over tt"*" — > ^jl^Vh- Moreover, the spectral 
distribution of neutrinos produced in the decay processes and the probability that a pion 
decays in flight in neutrinos show large deviations from the standard results. 
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1 Introduction 

The OPERA collaboration recently reported a 6.2 a evidence for superluminal neutrino 
velocity [1]. The reported value S = v -1 = (2.37 ± 0.32 (stat.) toU (sys-)) x 10^^ was 
obtained by observing the arrival times of the muon neutrinos of the CNGS beam, that 
travel over a baseline L ~ 730 km between CERN and Gran Sasso and have an average 
energy of about 17 GeV. 

One possible interpretation of the OPERA anomaly is that neutrino propagation does 
not obey Einstein special relativity. This interpretation should be confronted with the 
existing experimental constraints. On one side, the OPERA claim is compatible with the 
previous result of MINOS, which reported a deviation 6 = (5.1ib2.9) xlO~^ [2] for an average 
muon neutrino energy around 3 GeV. On the other side, if neutrino dispersion relation does 
not depend on the neutrino flavor, as it is suggested by neutrino oscillations, and if it is 
independent of energy, the deviation reported by OPERA is not compatible with the bound 
6 < 2x 10^^ [3] which has been deduced on the basis of the SN1987A neutrino observations 
[4], where the relevant distance and neutrino energy are 51 kpc and 20 MeV, respectively. 
The SN1987A bound can be avoided by assuming that the neutrino dispersion relation has 
a non-trivial energy dependence. As an example, one can assume that 6 scales as a power 
law of the neutrino momentum 5 ~ [5~7] with a sufficiently high exponent; alternatively, 
one can postulate that 6 is zero at SN1987A energies and has a sharp transition to the 
value observed by OPERA in the region ~ 0.1 — 1 GeV [6]. Stringent constraints also 
arise from the process of electron-positron pair creation — t- i/^ e~^e~ that rules out most 
of the above proposals [8-11]. In addition, other processes and most notably the very 
well measured charged pion decay process are expected to yield significant constraints on 
the superluminal nature of neutrinos. Kinematical effects due to non-standard neutrino 
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propagation were considered in [12] and [9, 13]. A calculation of the pion decay rate in 
the pion rest frame was performed in [14] , where energy-independent modifications of the 
velocity were considered and it was assumed that Lorentz violating terms are the same for 
neutrinos and charged leptons. 

Charged pion decays are the basic process for high-energy neutrinos production, and 
have a prominent role in long baseline experiments, in the production processes of atmo- 
spheric neutrinos, etc. For this reason, we focus on the study of charged pion decays and we 
determine a general expression for the decay rate of charged pions, assuming non standard 
neutrino dispersion laws. We adopt a "bottom-up" approach: starting from the dispersion 
law of neutrinos and assuming that most of the basic symmetries are not broken, we derive 
the pion decay rates. In this approach we preserve rotational invariance and space-time 
translational invariance, but we assume that boost invariance is broken. We also assume 
that the principles of quantum mechanics are valid for neutrinos and that the standard 
weak interactions are the same as in the standard model. Then, we consider various neu- 
trino dispersion relations proposed in the literature in connection with the OPERA result 
and we show that in each of the different cases, the decay rates, the decay probability of 
pions as well as the spectral distribution of neutrinos change significantly as compared to 
the standard case. 

We do not attempt the construction of a complete theoretical model that encompasses 
the various experimental results. Such a model should be able to reproduce the SN1987A 
data as well as the results of the OPERA experiment in a consistent framework. Our 
goal is instead, to construct a simple effective framework that allows us to investigate and 
eventually to falsify the neutrino dispersion relations recently proposed by various authors 
and/or the basic assumptions adopted in our approach. 

This paper is organized as follows. In Sect. 2, we recall the connection between neutrino 
velocity and dispersion relation. In Sect. 3, we write the modified Dirac equation assuming 
a generic velocity dispersion, and various useful identities. In Sect. 4, we recall well- 
known facts on neutrino interaction that remain unchanged in the present context. Then, 
in Sect. 5, we evaluate the decay rate for the vr"*" — )• fi'^u^ and vr^ — J- e'^i'c processes for a 
general dispersion relation. In Sect. 6, we discuss the implication of the different dispersion 
relations on the spectral distribution of neutrinos, on the branching ratios of the different 
decay modes and on the total decay probability. Finally, we summarize our results in 
Sect. 7. 

2 The basic assumptions 

Let us assume that it exists at least one reference frame in which space and time translations 
and spatial rotations are exact symmetries. This implies that energy, momentum and 
angular momentum in that reference frame are conserved, i.e., we can write for a generic 
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physical process: 

Eini = -E^fin 

Pini = Pfin (2.1) 
Jini — Jfiii 5 

where E'ini (-Efin), Pini (Pfin) and jmi (ifin) are the total initial (final) energy, momentum 
and angular momentum in the process. 

In this particular reference frame, that we identify with the laboratory frame, the 
energy of a given particle cannot depend on space and time coordinates and it can be only 
a function of the modulus of its spatial momentum. The specific form of the dispersion 
relation may depend on the particle type. However, for all particles except neutrinos, 
there are very strong constraints on possible deviations form the standard dispersion law, 
as reviewed in [15]. Therefore, we assume that all particles except neutrinos satisfy the 
well-known relation between energy and momentum that is provided by special relativity: 



= + ^2 ^ (2.2) 

where m[ is the particle mass. 

To infer the dispersion relation of neutrinos, we consider that the velocity u of a given 
particle is related to its energy by the Hamilton-Jacobi equation 

^ dE p 

V = — - , (2.3) 
dp p 

where p = \p\- In wave mechanics, the above expression gives the group velocity of the 
particle wave packet. This means that if we know the neutrino velocity v{p) as a function 
of its momentum we can determine its energy by performing an integration: 

E{p) = r dqv{q) + EQ. (2.4) 
Jo 

The available experimental data yield useful phenomenological constraints: 

1. The observation of neutrino flavor oscillations puts very strong bounds on the possi- 
bility that neutrino dispersion relations depend on their flavor [7]. Thus, we assume 
that the Lorentz violating effects are flavor universal. 

2. Assuming that neutrino and anti-neutrinos have the same velocities, we can use the 
bounds obtained from SN1987A to conclude that 

v(p) = 1 for p < 40 MeV, (2.5) 

with an accuracy at the level of 10^^ or more.^ This implies that at low energies we 
can safely neglect any Lorentz violating effect in our calculations, and set 

E{p)=p, (2.6) 



""^Ref. [3] quantifies j5| < 2 x 10^'' assuming that the neutrinos arrived within 3 hours from the time 
when the light arrived. This bound can be tightened by a factor of ~ 6 by modeling the propagation of the 
shock wave as in [16] and including in the analysis the observations of [17]. 
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where we set £Jo = in Eq. (2.4) and also neglect the neutrino mass. Indeed, stringent 
limits apply: \Eq\ is smaller than 5 eV from beta decay experiments [18], that imply 
also the bounds on neutrino masses of 2.3 eV [19] and 2.5 eV [20]; the analysis of 
SN1987A itself restricts the mass below 5.7 eV [21] (all bounds are given at 95% CL). 

3. OPERA findings suggest that neutrino velocity deviates from 1 at higher energies. 
Therefore, we write the velocity of neutrinos as 

v(p) = 1 + S{p) , (2.7) 

where 6{p) is an unknown function of momentum which characterizes the deviation 
from the standard value. According to the results of OPERA 

6{p) ~ ~ 2.5 X 10"^ for p ~ 17 GeV , (2.8) 

meaning that at some energy E{p) ^ p. 

For later convenience we also define 

F{p) = ^ = 1 + [\z 6{p z) , (2.9) 
P Jo 

where z is an integration variable. The function F{p) is the phase velocity of neutrinos and 
a measure of the deviation of the dispersion law of neutrinos from the standard expression; 
the deviation being proportional to the average of the function 5{p) over the interval [0,p]. 

In the rest of the paper, we discuss the implications of a generic F{p) ^ 1 on the 
neutrino wave function, on the description of neutrino interactions and on the pion decay 
processes. Then, in Section 6 we consider various specific expressions of F{p), motivated by 
the findings of OPERA collaboration, and determine the corresponding pion decay rates. 

3 The neutrino wave function 

According to the principles of quantum mechanics, energy and momentum are represented 
by the operators E — )• idt and p — t- — iV, respectively. The wave function of a neutrino 
with energy E and momentum p is thus described by 

i/j{x]p,X) = u{p,X)exp[—iPu-x] (3-1) 

where x = {t,x) and 

Pu^{E,p) . 

The spinors u{p, A) describe the intrinsic properties of the wave function and when E = p, 
i.e., F{p) = 1, they obey the standard Dirac equation. When F{p) ^ 1, the Dirac equation 
should be modified as follows, 

[7°E-F(p)7-p] n(#,A) =0 (3.2) 
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where = (7*^,7) are the usual Dirac matrices. Indeed, by using the anticommutation 
properties of the 7^^, we immediately find that the positive energy solutions of the above 
equation satisfy: 

E{p) = F{p)p, (3.3) 

as it is required by Eq. (2.9). Note that, even if Eq. (3.2) is multiplied by a non-zero 
multiplicative factor, e.g., its solutions and Eq. (3.3) remain unchanged. 

For the generic neutrino wave function ip{x), the equation of propagation reads, 

7°at + F(|V|)7-v] V(2;) = 0, (3.4) 

where |Vp = V • V. This is equivalent to Eq. (3.2) for plane waves and follows from the 
Lagrangian density, 

C = ii){x) [7°^* + F(|V|) 7- V] i^{x), (3.5) 

The breaking of SU(2) gauge symmetry and dependency of F on momentum will be dis- 
cussed in the next and subsequent sections respectively. 
For notational convenience, we define the quantity 

Pv = ip,P) (3.6) 

that coincides with the energy-momentum vector of a massless particle with momentum p 
in special relativity. By using this definition and considering that E = F{p) p, the condition 
in Eq. (3.2) can be rewritten as: 

u{p, A) = (3.7) 

where we adopted the usual convention ^ = 7^a^. This condition coincides with that 
derived in special relativity for a spin 1/2 particle with vanishing mass, indicating that 
the rotational properties of neutrino wave function in Eq. (3.1) are identical to those of 
standard neutrinos. It is convenient to require that the spinors u(p, A) have a definite 
helicity, i.e. 

u(p, A) = A u(p, A) with A = ±1 (3.8) 

p 

and we normalize the spinors in such a way that 

u{p,X)ru{p,X) = 2{p,r 
which implies that the density matrix is given by 

A 

4 Neutrino interactions 

In the standard model, weak interactions are due to the coupling of quarks and leptons 
to W and Z bosons, described by the charged-current and neutral-current interaction 
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Hamiltonian density 

^F"^ = :^j^''w^ + h.c. (4.1) 

Z COS t'w 

where g is the SU(2)i couphng constant, is the weak angle and the charged and neutral 
weak currents are given by: 



■cc 

J 1^1 



2 ^ 7m + • • • (4-3) 

£=e, fi, T 

j^^ = ^ 7/. i^tt + ■ • • (4.4) 

£=e, fi, T 

where i^iLix) = [(1 — 7^)/2] f£(x) and £l{x) = [(1 — 7^)/2] £{x) are the left-handed neutrino 
and charged lepton fields respectively, and we have written explicitly only the terms involv- 
ing neutrinos. Notice that, although we assume that the interaction vertices of neutrinos 
are as in the standard model, given the non-standard dispersion relation of neutrinos the 
space and time dependence of the interaction Hamiltonians necessarily change. 
The neutrino free field operators in the interaction representation are given by 



(27r)32p 



bi{p,-l)i^ix;p,-l) + 4ip,+l)i^^{x;p,+l)\ , (4.5) 



where the functions ij){x]p, A) have been defined in Eq. (3.1). The charge conjugate function 
Tp^ = C'lp* (i.e., the positive frequency component) is commonly rewritten introducing an 
auxiliary spinor f (p, A) as follows: 

^p^{x;p, A) = r]{p, A) v{p, A) exp [+iPu ■ x] , (4.6) 

where r/ is a conventional phase. From this relation one can check the properties of the 
auxiliary spinor, obtaining in particular that its density matrix is the same as in Eq. (3.9): 

v{p, X)v{p, A) = p^. The operator bi{p, —1) {dl{p, +1)) in Eq. (4.5) destroys (creates) a 
neutrino (antineutrino) of flavor i with momentum p in a negative (positive) helicity state. 
These operators are normalized to give (0| vl{x) \p) = tp{x;p, —1) for the transition from 
the one-neutrino state \p) = h\p) |0) to the vacuum.^ 

Note that we are proposing a modification to the standard model Lagrangian in which 
only the neutrino propagation is changed, to account for the finding of OPERA. Henceforth, 
we shall refer to Eq. (3.5) as the minimal modification of the standard model Lagrangian. 
In principle, one could be interested to consider different modifications of the weak inter- 
actions, when the interaction vertices are also modified [22] or when the charged leptons 
propagate in non-standard fashion [14]. These schemes also have an impact on pion de- 
cay [14, 22], but they go beyond the minimal modification that is needed to account for 



In other words, we have the anticommutation relations 
similarly for the antiparticles. 
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OPERA findings; we shall not elaborate further in this direction. Finally, it is to be taken 
into consideration that the minimal modification of the standard model Lagrangian, which 
we propose does not correspond to the principle of 'minimal coupling', that would imply 
that the free Lagrangian of the left electron is the same as the one of the neutrino in the 
massless limit. 

Incidentally, note that the effective Lagrangian of Eq. (3.5) can be derived from an 
SU(2) symmetric lagrangian after spontaneous symmetry breaking. It is sufficient to endow 
the standard Lagrangian density with the gauge invariant term, 

6C = i *£(F-1)7-V*£ (4.7) 

£=e, fi,T 

where the auxiliary SU(2) singlet field is defined as. 




when (H^) 7^ 0, we have = f£, reproducing Eq. (3.5). Note that the covariant derivative 
coincides with the ordinary derivative in Eq. (4.7), and the interactions stay unchanged. 
Alternatively, one could modify the temporal part i^^^^d^ (or both the temporal and 
the spatial parts), still obtaining the same neutrino wavefunctions and thus deriving the 
same consequences for pion decay. For an estimate of the radiative corrections and a 
discussion of a plausible impact of this assumption on the charged leptons, see [7]. 

5 Charged pions decays 

As it was suggested in [9, 12, 13], the charged pion decay processes are extremely sensitive 
to possible modifications of the neutrino dispersion law. At sufficiently low energies, these 
processes are described by the usual effective Hamiltonian density: 

^efi = ^ (5V^) i;:£pt + h.c. With ^^^cosOcf. (5.1) 

where (p-,^ is the pion field, j^f^p^ represents the leptonic part of the charged weak current, 
Gf is the Fermi constant, Oq is the Cabibbo angle and is the pion decay constant. 

We calculate the differential decay rate of the process vr^ — t- i^i'e in the laboratory 
frame using: 

dr = ^|A^P, (5.2) 

where E^^ is the pion energy, Ai is the matrix element (i.e., the amplitude) of the considered 
process and the phase space factor is defined as: 
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Note that the phase space differs shghtly from the conventional expression. Indeed, 
the standard volume element d?p/i2E) is replaced in our Eq. (5.3) by the factor d^p/{2p), 
which follows from the normalization of the neutrino spinors given by Eq. (3.9).'^ 
By using the neutrino field operator, Eq. (4.5), we obtain: 

\M\^ = eT^[M^-l5){^i + me)f'M-l5)p,] , (5.4) 

where Pu is defined in Eq. (3.6) and coincides with the energy-momentum vector of a 
massless particle in special relativity. The invariant amplitude in Eq. (5.4) is formally 
identical to that calculated for a massless neutrino in the standard theory, however when 
the energy and momentum conservation laws are used, the invariant amplitude differs 
from the standard expression. Indeed, in the phase space element, Eq. (5.3), the neutrino 
momentum P^, = (E, p) and not pu = {p, p) does appear. 

As a matter of fact, using the standard spinor algebra and the condition Pi = P-^ — Py 
from energy-momentum conservation, we obtain: 

|7W|2 = se [{mj - ml) P^ ■ p, + ml P, ■ p,] (5.5) 

where we introduced the definition 

ml{p) = E"^ - p^ . (5.6) 

We, then, can use the following relations: 

P^-p^ = p{E-p) (5.7) 
P.-Pv = \ [ml-mj + {E-p){E+p-2E^)\ (5.8) 

to cast the result into the form: 

\M\'' = Afimj ml + {E- pfpl - [mj - [E - p){E + p - E^)f] . (5.9) 

Setting E = p m. the above, we recover the well-known standard result 4^2m|(m^ — mf), 
that vanishes in the chiral limit m^ — )■ 0. Note that the same happens to the first order 
term in {E — p), namely 8^'^{E — p)m'j{E + p — E^^). 

The chiral limit remains important also in the non-standard case E ^ p, and it is 
particularly useful to understand the main features of the process vr"*" — )• e+z^e- Indeed, 
note that: 

1. If the pion is at rest [i.e., pn = 0), momentum conservation implies pi = — and 
setting m£ — )• 0, we have E-,^ = E + p and then \M\'^ — t- 0. 

2. The above does not hold anymore as soon as p-,^ ^ 0. When the pion decays in 
motion, we can rewrite |A^P = S^'^{E — p)'^{pep + Pi ■ p), that is evidently positive 
and second order in {E — p). 



^ Since the phase space is multiphed by the probabihty and the amphtude of emission M is 

proportional to the wavefunction of the emitted particle, the adopted normalization convention does not 
affect the final result, as it should be for any consistent calculation of the rate. 
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This has the important consequence that the contribution of Tr~^ — )• e^Ve to the total decay 
rate of moving pions can become large (even dominant) in non standard scenarios. 

Upon substituting Eq. (5.9) in Eq. (5.2) and after some algebra, the total pion decay 
rate can be expressed as an integral over the neutrino momentum: 

dp e = e,n (5.10 
dp 

where 

^ = aJe^p^ i^"^' - - - + ml - EAmj - ml)]} (5.11) 

gives the spectral distribution of neutrinos in the final state. The integration limits Pmin 

and 

Pmax are determined by studying the kinematic of the decay process and can be obtained 
by solving the equations: 

'^1 + "T-^(Pmax) -mj = 2£'(pmax)-E',r " "^PnPiaax (5.12) 

ml + mlipmin) -mj = 2E{p^i^)E^ + 2p^pmin (5.13) 

where E{p) = F{p) p and = p^ {F{p)'^ ~ !)• I^i the standard case, one obtains the values 
Po,max = (1 - mj/ml) {Et, + p,r)/2 and po.min = (1 - my ml) [E^, - p-„)/2 that are shown 
as black dashed lines in Fig. 1. 

If we consider a generic neutrino dispersion relation, the above equations have to be 
solved numerically. However, it is possible to obtain an analytical expression for Pmax and 
Pmirn if we assume that the neutrino velocity is constant, i.e., y = 1 + 5. In this case, 

= (5 (2 + (5)p^ and (E — p) = 6p that allows us to rewrite the conditions (5.12), (5.13) 
as second order equations in the neutrino momentum, and the upper limit pmax is given 
by the relation: 

(l-^) x{E^+p^) 

Pmax = ; (5.14) 



while Pmin can be obtained replacing + — )• i?^ — = ml/(ET^ + P-k)- 

Let us analyze this expression in the chiral limit mi — )• 0, when it takes a very trans- 
parent form. If the pion momentum is sufficiently high and if 5 > 0, the upper integration 
bound changes as follows: 

d{E^+Pjr) ^y5{2 + 5) 

From Eq. (5.15), one sees that Pmax tends to zero for large values of p-,^. Thus, the phase 
space for pion decay is strongly reduced with respect to the standard case. In all other 
cases {i.e., small pion momentum, or 5 < 0) we have instead 

Pmax= , (5.16) 
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that is a just a minor modification of tfie standard expression. 

Two final remarks are in order: (i) In the case when 6 > 0, the condition on the pion 
momentum of Eq. (5.15) can also be expressed as 



rriuipTr) > m- 



(5.17) 



namely as a condition on the minimum 'effective' neutrino mass defined in Eq. (5.6). 
{ii) Similar reasoning applies when we consider deviations from the chiral limit, as one can 
verify from the full expression of Eq. (5.14). In particular, for (5 > 0, the region where the 
Lorentz violating effects become relevant can still be identified by Eq. (5.17). 

6 Applications 

The assumption that neutrino propagation deviates from the expectations of Einstein rel- 
ativity implies that the dependency of the lifetime and the branching ratios on the pion 
momentum differ from the standard expectations. In the following we quantify these de- 
viations, considering specific dispersion relations that have been proposed in the literature 
and comparing the results with those obtained with the standard dispersion law, E = p. 

Dispersion relations The neutrino dispersion relation that have been proposed in con- 
nection with the OPERA results are the following ones: 

A. The neutrino velocity exceeds the velocity of the light by a constant multiplicative 
factor. In this assumption 



with = 2.5 X 10 ^ to match OPERA results, but in disagreement with the SN1987A 
constraints. 

B. The behaviour of neutrino velocity can be parameterized by a power law, i.e., v(p) — 
1 oc p". In this assumption, we have: 



where a > 3 to avoid the SN1987A bounds at lower energy. We take a = 4 in our 



C. The neutrino velocity behaves as a step transition from v = 1 at low energies, to the 
value measured by OPERA v = 1 + (5, at a transition momentum pt ~ 100 MeV. 
This corresponds to assuming: 



F{p) -1 = 5 






analysis. 




when pij < pt GeV 
otherwise 
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Log||,(;VlGeV) 



Log,„(p,/lGeV) 



Figure 1. The kinetic limits p-^^^y^/ E„ andpyj^^^/ for the processes 7r+ — > p^v^ and — > e^v^ as 
a function of the pion momentum p^r- The black dashed lines correspond to the standard case. The 
coloured lines correspond to the values of p^^^/ E-j^ obtained with neutrino dispersion laws discussed 
in Sect. 6. The effects of the considered modifications of the lower limit p^-m/ Ett are not appreciable 
in this plot. 



D. Finally, we consider the extreme proposal from [10] which was built ad hoc to explain 
OPERA result while suppressing the neutrino pair emission process — )• e"*"e~. 
This corresponds to assuming: 



v(p) 



1 



where 5 



2.75 X 10-3 and 62 

5^ (1_ O.lGeV^ 



5i when p = 0.1 — 1 GeV 
62 when p = 1 - 100 GeV 
otherwise 

2.5 X 10^^, that leads to 



(6.4) 



F{p) -1={ 



62 




1 



p 

IGoV 



+ 61 



IGcV 



0.1 GcV 

P 



when pi, = 0.1 — 1 GeV 
when p^ = l- 100 GeV (6-5) 
otherwise 



In this case the dispersion law is intentionally constructed to have F{p) < 1 for any 
value of p. The effect on F{p) produced by the positive 5 = 2.5 x 10~^ required by 
OPERA data has been cancelled by artificially postulating that 6 
the region p = 0.1 — 1 GeV. 



-2.75 X 10-3 in 



Next, we discuss the effect of these dispersion relations on pion decay kinematics, spectral 
distribution of muon and electron neutrinos, as well as the effect on the pion decay lifetime. 

Kinematics In Fig. 1 we show the kinematic limits for the two channels vr^ — t- h'^v^ and 
vr^ —7- e'^i'e- The dashed lines correspond to the standard case. The solid coloured lines 
represent the values of Pmax/E-,^ and Pmin/E-,^ obtained with the above mentioned neutrino 
dispersion laws. We see that in all cases the kinematics of the process is radically affected 
at energies E^^ > 10 GeV. For the dispersion laws corresponding to the cases A, B and C, 
this result can be understood by noting that the non-standard terms in Eq. (5.12) scale 
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Figure 2. Left panel: The spectral distribution of muon neutrinos produced in the process 7r+ — > 
fJ-^i^^i by pions with momentum = 40 GeV. Right panel: The spectral distribution of electron 
neutrinos produced in the process 7r+ C^Vg by pions with momentum p-^ = 40 GeV. The labels 
A, B, C and D refer to the scenarios described in Sect. 6 . 

approximatively as ~ 25 and (E — p) ~ 6 p. If we consider that 6 ~ 2.5 x 10~^, as 
indicated by the OPERA experiment, and we use pmax ~ Po.max as a rough estimate, we 
can calculate that non standard terms become dominant when ~ rxiT^l^fb ~ 20 GeV. In 
the case A, it was already noted by [13] that the modification of the neutrino dispersion 
relation produces a reduction of the phase space for pion decay. This is in agreement with 
our results and analytical expressions, see Eq. (5.14) and the subsequent discussion. We 
find a similar behaviour in the case B and C where we have -F(p) > 1 and rhKp) > 0. 

In the case D, the dispersion law leads to a more complicate behaviour. As shown in 
the right panel of Fig. 1, for the process tt~^ — )• e"*"fe there is no appreciable effect on the 
integration bounds. However, for the process 7r+ — )• /u+f^ the phase space available for the 
decay is much larger in the interval 1 GeV < p < 100 GeV. In particular, for p ~ 100 GeV 
we find that pmax can become as large as 0.9 Et^. Notice that the dispersion law of the case 
D was built to take into account the results of the OPERA experiment and to suppress 
the pair creation process t'^ — t- e~^e~. However, it produces a big change on the phase 
space of the vr"*" — t- fi'^v^ process and, as we shall see, strongly modifies the decay rates. 

Momentum distribution In Fig. 2, we show the momentum distribution of neutrinos 
(normalized to one) produced by the decay of pions with momentum p^^ ~ 40 GeV. This 
was calculated according to our description of the decay modes vr"'" — )• ^^f^ (left panel) and 
TT^ — )• e'^i'e (right panel) given in Eq. (5.11), or equivalent, in Eq. (5.9). The predictions 
obtained in the various scenarios are completely different among each others. The very pe- 
culiar spectral shapes for vr^ — )• e^fg can be understood considering that Lorentz violating 
terms represent the dominant contribution to the decay processes, while in the standard 
scenario - full line in Fig. 2 - the decay process is suppressed by chirality arguments. For 
both decay modes, the differences w.r.t. the standard expectations are large. This shows 
that the measurement of spectral distribution of neutrinos (or, equivalently, muons and 
electrons) produced in pion decay is a sensitive tool to probe and possibly to falsify the 
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Figure 3. Modifications of the decay rate of pions and of the pion lifetime with respect to the 
standard case. The labels A, B, C and D refer to the scenarios described in Sect. 6. Left-upper panel: 
The ratio T^/T^ between the rate o/tt"*" — > pL^v^ obtained with modified neutrino dispersion laws and 
that obtained with the standard dispersion law. Right-upper panel: The ratio Tc/Tc between the rate 
o/ 7r+ — >■ e+i-'o obtained with modified neutrino dispersion laws and that obtained with the standard 
dispersion law. Left-lower panel: The ratio between the decay rates of the processes tt"*" — >■ e'^Vc and 
7r+ — > fJ.^h'fj,. Right-lower panel: The pion lifetime as function of the pion momentum. Note the 
logarithmic scale of the plots. 

assumed neutrino dispersion law and/or the set of hypotheses adopted in our approach. 

Decay rates and lifetime In the upper panels of Fig. 3, we show the ratios F^/F^ and 
Fe/Fe as a function of the pion momentum , where F^^g are the decay rates of vr^ — t- fi^v^, 
TT^ — 7- e'^Ve, and F^^e represents the standard model predictions. For vanishing values of 
Ptt, the effect of Lorentz breaking terms is proportional to 5, in agreement with [14], and 
is not appreciable in Fig. 3. For larger values of p-j^ the effect is much larger, because the 
breaking of Lorentz invariance alters the normal scaling F^ cx 1 / producing very peculiar 
behaviours with energy. In particular, for i?^ ~ 100 GeV, that we take as a rough estimate 
of average energy of pions produced in OPERA, the pion decay rate to muons is decreased 
by a factor ~ 1/50 in the scenarios A and C, and by a factor ~ 1/3 in the scenario B, while 
it is increased by a factor ~ 200 in the scenario D. Even more significant effects are obtained 
for the TT^ — )• e'^fe decay process. The rate of this process is enormously enhanced due to 
the Lorentz violating terms in the matrix element. For a moving pion (i.e., Ptt 7^ 0), the 
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Figure 4. Left Panel: The probability that a pion produces a muon neutrino in a tunnel of length 
L = 1 km as function of the pion momentum. Right Panel: The probability that a pion produces an 
electron neutrino in a tunnel of length L = 1 km as function of the pion momentum. The labels A, 
B, C and D refer to the scenarios described in Sect. 6. 

matrix element and hence the decay rate does not vanish in the chiral hmit (see Eq. (5.9) 
and related discussion), with the important consequence that the process vr^ — ?■ e'^v^ can 
provide a non negligible contribution to the total decay rate. This can be seen in the left- 
lower panel of Fig. 3 where we show the ratio Tg/r^ as function of the pion momentum. 
In the cases A and C, we obtain Fe > F^ in the interval = 20 — 40 GeV, much larger 
than the experimentally value observed in the pion rest frame Fg/F^ = 1.2 x 10~^ [23-25]. 
In the case D, we have Fg ~ F^ for E^^r > 10 GeV. The change in the matrix element, 
Eq. (5.9), as well as the change in the pion decay kinematics, portrayed in Fig. 1, clearly 
affect the total decay width, and hence the lifetime of pions that decay in motion: 



As shown in the right-lower panel of Fig. 3, this scenario implies significantly large devia- 
tions from the standard scaling of the pion lifetime, f oc E-^. In this respect, it is interesting 
to note that the old measurement of [26] found that the lifetime of the pion that decays in 
flight at = 300 MeV agrees with the time dilatation predicted by Einstein relativity at 
the ~ 0.15% level of accuracy. At this energy, the modifications expected in the cases A, 
B and C are small, whereas in the case D one obtains a ~ 0.5% reduction of the lifetime. 

Probability of decay The above results have important implications for OPERA and, 
more in general, for high-energy neutrino experiments, since the expected neutrino signal is 
radically changed both in spectrum and in composition. In Fig. 4 we show the probability 



that a pion of a given momentum produces a neutrino ui in a tunnel of length L = 1 km 
(note that the standard and the additional dependences of the lifetime and of the widths on 
the momentum are fully included). The probability is suppressed for momenta above 



1 



(6.6) 



Fo(p7r) + ^fiiPw) 




(6.7) 
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Ptt ~ 20GeV in the scenarios A, B and C. Instead for the case D, the suppression starts 
at much lower energy, due to the competing electron neutrino production. For the case A, 
the probability of decaying into an electron Ve (right panel of Fig. 4) is larger than the 
probability of decaying into a muon in the interval p-,^ ~ 20 — 40 GeV. In all the cases 
considered, the electron neutrino fraction in the beam is drastically increased with respect 
to the standard expectations, showing that the electron-to-muon decay rate, in the range 
relevant for the OPERA experiments, could be a sensitive probe for non standard neutrino 
propagation. 

7 Conclusions 

The recent OPERA experimental results demand for non standard neutrino propagation. 
In this work, we discuss some implications of this assumption. We show, in particular, 
that within a set of well-defined hypotheses (i.e., the space and time translational invari- 
ance, the rotational invariance, the basic quantum mechanical principles and the standard 
weak interaction hamiltonian) it is possible to calculate the consequences of non standard 
neutrino propagation on a generic physical process which involves neutrinos. 

We then apply our approach to the charged pion decay processes vr"'" — )• ^'^v^ and 
vr"*" — )• e^fg. We consider various neutrino dispersion relations which have been proposed 
in connection with the OPERA result. Namely, we assume that: the neutrino velocity 
is constant (case A in the text); the neutrino dispersion law scales as power-law of the 
neutrino momentum (case B); the neutrino velocity has a sharp transition at ~ 100 MeV 
to the OPERA measured value (case C); the dispersion law is chosen ad hoc in order 
to suppress the neutrino pair production process (case D). The impact of the assumed 
dispersion laws on the decay kinematics is shown in Fig. 1; the modification of the spectral 
distributions are shown in Fig. 2; the modification of the decay rates and of the pion lifetime 
are given in Fig. 3; finally, the impact of the various dispersion laws on the probability to 
produce electron and muon neutrinos in OPERA experimental setup is shown in Fig. 4. 

We conclude that, for all of the considered dispersion relations, the pion decay processes 
suffer a drastic departure with respect to the standard scenario in the energy interval 
relevant for OPERA. To quote a few eloquent numbers, the rate of vr^ — t- ix'^i'^ at E^^ = 100 
GeV decreases by about 1/50 and 1/3 in the cases A, C and in the case B respectively, 
while it increases by a factor ~ 200 in the case D. Moreover, in all the considered cases, the 
probability to produce electron neutrinos at the energies relevant for OPERA is drastically 
increased with respect to the standard expectations. 

All this shows that the study of the charged pion decay can be used as a sensitive 
probe to investigate and possibly falsify the assumed neutrino dispersion laws and/or the 
basic assumptions adopted in our approach. 
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